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SUMMARY 

A  technique  is  developed  for  calculating  rocket  base 
heating  and  spacecraft  heating  environments  due  to  particle 
radiation  from  a  single  nozzle  rocket  exhaust  plume.  The 
technique  has  proved  successful  when  applied  to  a  single 
nozzle  exhausting  into  a  rarefied  atmosphere  on  the  basis  of 
comparison  of  predictions  with  experimental  results. 

The  analysis  treats  radiation  from  a  cloud  of  particles 
as  that  from  an  equivalent  radiating  surface.  Thus,  the 
problem  is  reduced  to  the  determination  of  the  proper  values 
of  the  apparent  surface  e  miss  ivity  and  the  effective  tempera¬ 
ture.  The  particle  flow  field  information  (particle  concen¬ 
trations,  temperatures,  and  trajectories)  necessary  to  de¬ 
termine  these  two  quantities  is  provided  by  a  two-phase  flow 
field  computer  program  developed  by  STL. 

In  defining  the  apparent  emissivity  of  the  particle  plume, 
an  analogy  with  neutron  scattering  for  a  cylindrical  cloud  is 
adopted  which  shows  the  apparent  emissivity  to  be  dependent 
on  particle  emissivity  and  cloud  optical  thickness.  Since  the 
plume  is  non-uniform  in  particle  size,  concentration,  and 
temperature,  certain  averaging  techniques  are  used  to  define 
mean  values  of  optical  thickness  and  temperature. 
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NOMENCLATURE 


Symbols 

2trr 

a  = 


E 

•s / 

E 

F 


I 

k 

L 


i 

m 

N 

P  i  - 1 ,  k 
Pn 

Pc 

Pe 

Q 

q 

r 

T 

x 


z 

a 

p 

V 


2 

photon  (radiant  energy)  flux,  Btu/ft  -sec 

integroexponential  function  (Equation  18) 

shape  factor  from  unshaded  portion  of  the  particle  plume 
to  a  surface  element 

radiant  intensity,  Btu/ft^-sec-steradian-micron 

cosine  distribution  exponent  (Equation  7) 

total  length  of  path  across  a  particle  plume,  ft 

distance  along  view  path  from  a  surface  element  to  the 
particle  plume  boundary,  ft 

length  through  a  particle  plume,  ft 

number  of  plume  segments 

3 

particle  concentration,  particles /ft 
collision  probability  (Equation  2) 

^n,  k 

rocket  chamber  pressure,  psia 
ambient  pressure,  psia 
efficiency  factor  (Equation  23) 
radiant  heat  rate,  Btu/ft^-sec 
particle  radius,  ft  or  microns 
absolute  temperature,  °R 

distance  along  path  from  the  boundary  of  a  particle  cloud,  ft 

x  cos  0  (Figure  3) 

absorptivity 

neutron  blackness  (Equation  2) 

angle  of  the  plane  of  view  with  respect  to  the  plume  axis, 
degrees 
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NOMENCLATURE  (Continued) 


Symbols 

c  =  emissivity 

t)  =  temperature  averaging  factor 
0  =  directional  angle,  degrees  (Figure  3) 

X  =  wavelength,  ft  or  microns 
fx  =  cos  0 

7T  =  infinite  product  symbol 

p  =  reflectivity 

<r  =  cross  section,  ft^ 

2  o  4 

<r,  =  Stefan-Boltzman  constant,  Btu/ft  -sec  -  R 

b 

r  =  optical  thickness 

</>  s  dummy  variable  in  integroexponential  function; 
represents  t  as  used  herein 

Subscripts 

a  =  apparent 
i  =  dummy  number 
m  =  number  of  plume  segments 
M  =  maximum 
n  =  dummy  number 
o  =  initial  or  incident 
p  =  particle 
s  =  scattering 
t  =  total  or  extinction 
tr  =  transmitted 
a  =  absorbed 
X  a  function  of  wavelength 

Superscripts 

—  a  average  or  effective 
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I.  INTRODUCTION 


The  employment  of  single  nozzle,  aluminized  composite  propellant 
rockets  to  boost  suborbital  payloads  or  to  inject  spacecraft  into  final  orbits 
results  in  primarily  radiant  heating  from  the  rocket  exhaust  to  the  boost 
rocket  base  and  to  the  spacecraft.  It  is  important  that  the  heat  rates  be 
predictable  to  permit  efficient  design  of  thermal  insulation  and  shielding 
and  to  permit  evaluation  of  the  effects  of  such  heating  on  spacecraft  solar 
cells,  sensors,  etc. 


If  the  altitude  range  is  restricted  to  that  experienced  by  the  latter 
stages  of  boosters  and  above,  it  can  be  shown  that  heating  of  surfaces 
located  forward  of  the  rocket  nozzle  exit  due  to  exhaust  gas  radiation  is 
negligible  compared  to  heating  from  aluminum  oxide  particle  radiation. 
This  is  true  because  of  the  relatively  low  temperature  and  emissivity  of 
the  expanded  gases.  The  problem,  then,  is  reduced  to  that  of  predicting 
the  thermal  radiation  from  the  exhausted  particle  stream  or  “plume”  to 
the  surface  in  question.  A  technique  sufficiently  accurate  for  engineering 
design  purposes  has  been  developed  and  is  the  subject  of  this  paper. 

The  approach  has  been  to  treat  the  cloud  of  particles  which  forms  the 
plume  as  an  equivalent  surface  radiator  whose  shape  is  that  of  the  plume. 
To  account  for  variations  in  radiation  intensity  along  the  surface  of  the 
plume  as  seen  by  the  surface  element  being  heated,  the  plume  is  divided 
into  a  number  of  segments  by  planes  emanating  from  that  element  (Figure  1). 
The  increment  of  radiation  incident  to  the  surface  element  from  each  seg¬ 
ment  is  determined  by  the  familiar  relation 


Aq.  =  AF.  €  o',  T. 

n  i  a.  b  i 


(i) 


A  summation  of  radiation  from  all  segments  visible  to  the  surface 
element  yields  the  value  of  the  total  radiation. 

The  shape  factor,  AF.,  of  each  segment  to  the  surface  element  is  a 

1  \ 

function  of  the  surface  shape  of  the  particle  plume.  Morizumi  has  devel¬ 
oped  an  analytical  method  for  calculating  shape  factors  to  cones,  hyper¬ 
boloids,  and  segments  thereof,  which  has  been  employed  to  arrive  at  the 
values  used  herein. 
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Figure  1.  Sketch  of  Nozzle  with  Gaseous  and  Particle  Plumes,  and 
Surface  Element  of  Base  Region 


The  remaining  two  quantities  to  be  determined  in  Equation  (1),  c  and 
T.  ,  are  functions  of  not  only  the  surface  shape  of  the  plume,  but  also  the 
aluminum  oxide  particle  concentrations  and  temperatures.  The  methods 
used  to  determine  these  two  quantities  are  discussed  here  in  detail. 

Techniques  used  to  arrive  at  all  three  quantities  rely  upon  particle  flow 
field-information  (particle  trajectories,  concentrations,  and  temperatures). 
The  flow  field  is  described  by  a  two  phase  flow  computer  program  developed 
by  Kliegel  and  Nickerson.^ 
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II.  ANALYTICAL  METHOD 


A.  APPARENT  EMISSIVITY 

1 .  Analogy  with  Neutron  Scattering 

In  determining  the  apparent  emissivity  of  a  particle  cloud,  an 

analogy  with  neutron  blackness  (probability  of  neutron  absorption  by  a 

medium)  is  used.  The  neutron  blackness  utilizing  the  probability  approach 

3  4  5 

is  expressed  by  Stewart  *  and  Anthony  in  their  studies  of  the  neutron 
scattering  problem  as 


P  = 


TT  k(T) 

i=i  1  1»K 


(2) 


The  collision  probability  P^_j  ^  is  defined  as  the  probability  that  a  neutron 
which  has  experienced  (i-  1)  collisions  with  other  particles  in  the  medium 
will  experience  at  least  one  more  collision  therein.  As  shown  in  Refer¬ 
ences  3,  4,  and  5,  the  collision  probability  increases  monotonically  with 
optical  thickness;  the  latter  is  a  measure  of  capture  probability  and  is  de¬ 
fined  in  the  following  section. 

A  similar  probability  function  can  be  used  to  define  the  radiant 
energy  absorbed  by  a  cloud  of  particles.  In  this  case,  the  photons  in  the 
particle  cloud  are  considered  to  replace  the  neutrons.  As  illustrated  in 
Figure  2,  the  radiant  energy  absorbed  in  the  first  collision  of  photons  with 
particles  is 


Likewise,  the  energy  absorbed  on  each  successive  collision  is 


AE  =  E  P  PlP  (1  -  p  ) 

a  o  o  lKp  rp 


AE  =  E  P  P .P,p  4  (1  -  p  ) 
a-  o  o  1  2Kp  '  Kp' 


AE 


=  E  P  P. 
an  °  °  1 


Pn-1  C  »  '  <>.'• 
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Figure  2.  Probability  of  Photon  Absorption  in  a  Particle  Cloud 


The  total  energy  absorbed  by  particles  through  collisions  is  then 


E«  =  i?i  iE°i  ■  Eo<‘  -  Vfr  +  PoPl'p  +  •  •  •  +  PoPlP2  • '  •  Vp" 
♦  •  •  Eo<‘  -  "p> £  >p'  TT  Pi-1 

K  J=1  ^  1=1 


The  absorptivity  or  emissivity  of  the  particle  cloud  is  the  ratio  of  the  total 
absorbed  energy  to  the  total  incident  energy,  which  is  defined  as 


a  = 
a 


(E  \  oo 


j-1 


J 

TTP; 

i=l 


i-1 


(3) 


Equation  (3)  can  be  simplified  by  assuming  that  the  particles  are  opaque, 

i.  e p  =(!-<)•  Thus 

P  P 


oo 


(4) 


A  comparison  of  Equation  (4)  with  (2)  shows  an  analogy  between  neutron 
blackness  and  apparent  emissivity. 

For  both  neutrons  and  photons,  the  change  of  intensity  due  to 

absorption  and  scattering  is  a  function  of  particle  total  cross  section  and 

particle  concentration.  The  total  cross  section  is  a  measure  of  the  energy 

taken  away  by  an  individual  particle  and  consists  of  absorbing  and  scattering 

cross  sections  (a.  =  <r  +(r  ).  It  is  seen  from  Equations  (2)  and  (4)  that  the  ratios 
t  s  a 

of  the  neutron  absorption  and  scattering  cross  sections  to  the  total  cross 

section,^  /<r  )and(<r  /o\l  are  equivalent  to  the  absorptivity  and  reflectivity, 
at  st* 

c^andl  respectively,  of  a  particle  in  the  radiating  medium. 

The  above  expression  can  be  applied  to  any  cloud  geometry.  The 
effects  of  the  cloud  shape,  the  total  cross  section,  and  the  particle  con¬ 
centration  are  reflected  in  each  of  the  probability  functions  P.  ^ 

2.  Probability  Functions 

The  determination  of  the  first  few  probability  functions  is 
shown  below  for  an  infinite  slab  of  particles.  Consider  a  beam  of 
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rays  traversing  a  slab  in  the  direction  of  0  with  respect  to  a  line  per¬ 
pendicular  to  the  slab  (Figure  3).  Take  an  element  of  this  beam  of  length 
dx  and  cross  section  A  at  x  distance  from  the  boundary.  The  total  number 
of  particles  in  this  element  is  NA  dx  and  the  total  cross  section  of  these 
particles  is  NA  dx.  Then  the  fraction  of  energy  removed  in  the  element 
through  absorption  and  scattering  is 

=  No“t  dx 

The  rate  of  change  of  the  intensity  of  radiant  energy  may  be  determined  as 
follows 


I  + 


||lj  dx  =  I  -(N<rt  dx)  I 


-(Nat)  I 


Integrating  this  once 


-No*  x 

1  =  *o  e  ' 


(5) 


(6) 


where  Iq  is  intensity  at  x  =  0. 

If  the  usual  cosine  function  distribution  of  the  radiant  energy 
intensity  is  assumed  for  entrant  distribution,  then  the  initial  intensity  IQ 
may  be  defined  in  terms  of  a  maximum  intensity  as  a  function  of  the  cosine 
of  directional  angle  0. 


I  =  I  (cos  0)k  =  I  / 
o  m  m 


Substituting  this  expression  into  Equation  (6) 

k  -Nv 

I  =  Impk  e 


(7) 


Differentiating 


the  above  equation 


-No.1  u* 
c  m 


-Noyc 

e 


(8) 
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The  energy  loss  per  unit  length  along  a  path  may  be  defined  as 


Combining  this  equation  with  Equation  (8) 

,  -  k  -N<r  x 

jdE|  =  -2irImN(rthi  e  d(i  dx 

The  total  energy  loss  for  a  plane  wave  traversing  across  the  slab  of 
thickness,  t,  is  then  determined  by  integrating  the  above  differential 
equation  over  the  slab  thickness,  t,  and  the  solid  angle  Zrr  steradians  as 
follows 


The  total  energy  incident  to  the  slab  can  be  determined  as 

fl  i.  2ttI 

Eo  =  2.  Io  djx  =  2ulm  /  H  =  TfrJ^7 


(10) 


The  probability  of  photons  colliding  with  particles  on  their  first  flight  in 
the  slab  is 


P 


o,  k 


-N<r^  z  /p 
e  dp  dz 


(ID 


As  shown  by  the  above  integral,  the  probability  function  is  dependent 
on  the  particle  concentration,  total  cross  section,  and  the  cloud 
geometry. 

In  order  to  determine  the  second  collision  probability,  the 
probability  that  photons  reflected  isotropically  by  a  particle  at  a  distance, 
z,  escape  through  the  surface  z  =  t  without  suffering  further  collision 
must  be  defined  first.  (This  probability  is  denoted  by  P(z).)  If  the  in¬ 
tensity  of  the  radiant  energy  reflected  isotropically  at  %  =  z  is  denoted  by 
Ig,  then  the  intensity  of  the  energy  transmitted  from  z  =  z  to  the  surface 
z  =  t  in  the  direction  of  9  may  be  expressed  as 


9 


L  =  I  e 
tr  s 


-Nat(t  -  z )  /p 


The  total  energy  transmitted  from  z  =  z  to  z  =  t  is  then  determined  as 

-I  -No-  (t  -  z)/p 


E,  =  2ft  T  I  dp  =  ZttI 
tr  /  tr  r  s  f 

Jo  ^o 


r1  -wo-  v 

J  6 

Jo 


d\x 


(12) 


and  the  total  energy  scattered  at  z  =  z  toward  the  surface  z  =  t  is 


Thus 


E 

s 


I  dp  =  2ttI 
s  r  s 


(13) 


P(z) 


-Ncrt(t  -  z)/p 
e 


dfi 


(14) 


where  the  factor  of  1/2  accounts  for  the  half  space. 

Similarly,  the  probability  that  photons  reflected  isotropically  at 
z  =  z  escape  through  the  surface  z  =  0  can  be  defined  as 


P  (z) 


-No\  z  /p 
e  dp 


(15) 


Hence,  the  probability  that  once  scattered  photons  escape  through  the 
surfaces  z  =  0  and  z  =  t  is 


Po,k<‘-Pl,k>V‘'p<k 


r 


+  1)  Nff 


r  i  r1  k-i  -Natz/^ 

[P(a)+P'(*)J  I  p  e  dfx dz 


The  second  collision  probability,  therefore,  is 


P 


1  ,k 


=  1 


W-l  -Nff. 

p  e  dp  dz  ( 1 6) 


By  successive  generations  of  the  probability  functions,  the  third  collision 
probability  can  be  similarly  determined  as 
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fl  1  -N<r.(z -z')/n 

p  e  dp  dz '  dz  (17) 

>o 


The  equations  for  the  collision  probability  can  be  simplified  by  using  the 
integroexponential  functions,  the  values  of  which  are  readily  available. 
The  function  is 


E .<*)  =  f1  pi-“  e"v' **  dp 

-'O 


-2  -<p/ p 


(18) 


Thus,  Equations  (11),  (16),  and  (17)  become 


f  nv 


Po(k  =  (k  +  D  N<rtEk+1(N<rtz)  dz 


(19) 


Pl,k=  1  f  N<rtgk+1(NV>  p2(Natz)  +E2(NV-Natz)]  dz  (20) 

o,k  Jo 

'  (“r)  h  :  f  *  Ez(Nv  •  N<rtz,l 

\  /  o,k  1  ,k  Jo 

rt 

I  £^+j(No-t  z1)  fi^No^z  -  No^z')  dz'  dz  (21) 

Jo 


P2,  k  =  1 


It  is  apparent  from  the  above  equations  that  the  collision  probability  is  a 
function  of  (N<r^t)  and  k.  The  first  parameter  (N<r^t)  is  called  "optical 
thickness"  and  is  a  product  of  particle  concentration,  total  (extinction) 
cross  section  and  slab  thickness.  The  second  parameter  is  the  exponent 
of  the  cosine  function  which  represents  the  entrant  distribution  of  radiant 
energy.  If  an  isotropic  emission  is  assumed,  k  =  0  is  used. 
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3.  Variation  of  Apparent  Emissivity  with  Optical  Thickness 

and  Particle  Emissivity _ 

Using  the  foregoing  approach,  Anthony  and  Stuart,  respectively, 
developed  the  neutron  blackness  curves  for  the  infinite  slab  and  cylinder 
in  connection  with  their  neutron  diffusion  studies*  The  optical  equivalent 
of  these  curves  is  presented  in  Figure  4.  The  apparent  emissivity  for  a  con¬ 
stant  particle  emissivity  is  a  monotonically  increasing  function  of  optical 
thickness,  t.  For  large  values  of  t  (optically  thick  cloud),  the  emissivity  curves 
approach  asymptotically  a  limiting  value,  regardless  of  the  cloud  geometry. 

For  very  small  values  of  t,  however,  the  apparent  emissivity  differs  for  dif¬ 
ferent  geometry.  For  very  small  values  of  t,  a  linearization  of  the  apparent 
emissivity  function  results  in  the  following  approximations. 

c  =  lc  t  for  a  cylinder 

a  p 

€  =  Z(  r  for  a  slab 

a  p 

It  is  observed  that  the  emissivity  for  a  slab  increases  (with  t  and  c  ) 

D 

twice  as  fast  as  that  for  a  cylinder  in  the  region  of  very  small  t. 

Although  the  configuration  of  a  rocket  exhaust  particle  plume  re¬ 
sembles  a  cone,  the  emissivity  for  the  cylindrical  configuration  is  used 
here  because  the  values  of  the  probability  functions  for  a  cone  are  analyti¬ 
cally  intractable  and  because  the  axial  symmetry  of  the  cylinder  approxi¬ 
mates  reasonably  well  the  rocket  plume  geometry. 

4.  Optical  Thickness  and  Particle  Emissivity 

It  has  been  shown  that  the  apparent  emissivity  for  particle  clouds 
is  a  functionof  particle  emissivity  and  optical  thickness.  Therefore,  if 
these  quantities  are  known  for  a  given  cloud  of  radiating  particles,  then  the 
neutron  blackness  curves  developed  for  a  similar  geometry  may  be  used 
to  determine  the  apparent  emissivity  for  the  particle  cloud. 

a*  Optical  Thickness 

The  optical  thickness  for  a  cloud  of  particles  of  uniform  size 
and  uniform  concentration  is  defined  as 

t  =  N 
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°» 'AiiAissiwa  iNaavddv 


13 


Inasmuch  as  the  particle  plume  for  a  typical  aluminized  propellant  rocket 
consists  of  varying  size  particles  and,  furthermore,  the  particle  concen¬ 
tration  is  not  uniform  in  the  flow  field,  the  optical  thickness  cannot  be  de¬ 
termined  directly  by  use  of  the  above  expression.  In  computation  of  the 
particle  plume  flow  field  by  the  computer  program,^  five  distinct  particle 
sizes  were  used  to  simulate  the  actual  flow  field.  Consequently,  a  super¬ 
position  technique  which  sums  up  the  thicknesses  for  all  particle  sizes  has 
been  adopted.  The  equation  is 

T  =  I  "t  f  n  Nn(x)  dx  (22) 

n=l  n  Jo 


The  total  cross  section,  cr^,  of  a  particle  plume  may  be  ex¬ 
pressed  in  terms  of  the  particle  geometrical  cross  section  as 


or,  =  Q,  irr“ 

t  t  p 


(23) 


The  coefficient  is  called  "extinction  efficiency  factor"  and  is  discussed 
quite  extensively  in  References  6  and  7.  This  efficiency  factor  is  a  function 
of  the  particle  size  and  shape,  the  refractive  indices  of  the  materials  in¬ 
volved,  and  the  radiation  wavelength.  The  values  of  the  total  cross  section 

for  small  carbon  particles  determined  by  the  Mie  theory  are  presented  in 

* 

Reference  8.  For  sapphire,  values  are  given  in  Reference  9. 

According  to  the  Mie  theory,  the  extinction  efficiency  factor 
of  a  particle  can  be  expressed  as  a  function  of 


As  this  value  increases,  the  value  of  also  increases  for  small  values  of 
a.  If  the  imaginary  part  of  the  refractive  index  is  of  a  small  order  of 
magnitude  (small  absorption  coefficient),  oscillates  for  intermediate 
values  of  a.  However,  for  large  values  of  a,  the  oscillation  damps  out  and 
the  value  of  becomes  equal  to  2.  This  fact  is  shown  in  Reference  6  for 
numerous  materials. 


* 


Aluminum  oxide  of  single  or  few  crystals 
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The  value  of  a  for  the  particle  plumes  under  consideration  is 
greater  than  6  since  the  mean  particle  diameter  is  4  microns  and  the  wave¬ 
length  for  maximum  radiant  power  is  from  1.5  to  2.0  microns.  For  this 
value  of  a,  the  efficiency  factor  may  be  assumed  to  be  2.  Thus,  the  total 
cross  section  cr^  in  Equation  (22)  is  taken  as  twice  the  geometric  cross 
section  of  the  spherical  particles. 

b.  Particle  Emissivity 

As  can  be  seen  from  Figure  4,  the  particle  emissivity  is 
equally  as  important  as  t  in  determining  the  apparent  emissivity  of  a 
particle  cloud.  A  theoretical  approach  to  the  determination  of  particle 
emissivity  is  to  apply  the  Mie  theory  to  a  sphere  of  appropriate  size  and 
observe  the  values  of 


from  which  a  total  emissivity,  can  be  obtained.  The  values  of  Qq^ 
and  Qt^  given  in  Reference  9  for  aluminum  oxide  (sapphire)  result  in  a 
value  of  considerably  less  than  0.1.  Higher  values  of  particle  emis¬ 
sivity  determined  grossly  from  measurements  of  rocket  exhausts  are  in 
conflict  with  this  result.  It  is  suspected  that  the  optical  properties  for 
sapphire,  upon  which  the  calculations  of  Reference  9  were  based,  are  not 
representative  of  those  for  aluminum  oxide  particles  in  rocket  exhausts. 
One  argument  which  has  been  advanced  is  that  the  absorptive  properties 
of  the  particles  are  higher  than  those  of  sapphire,  due  to  a  polycrystalline- 
porous  structure  of  the  former.  Lee  and  Kingery^  have  explored  this 
argument  in  some  detail.  The  high  experimental  values  of  particle  emis¬ 
sivity  seem  to  justify  the  assumption  of  opaqueness  in  Equation  (4). 

For  application  to  the  present  problem,  a  single  value  of 
=  0.25  was  chosen.  This  value  was  determined  based  on  extrapolation 
of  bulk  aluminum  oxide  emissivity  measurements .  The  single  value  is  justi¬ 
fied  because  of  the  small  sensitivity  of  particle  emissivity  to  temperature 
in  the  range  of  the  temperature  of  the  particles.  The  values  presented 
and  discussed  in  Section  III  tend  to  confirm  the  choice  of  c  =0.25. 


15 


The  mechanics  used  to  arrive  at  the  particle  emissivity  from 
rocket  plume  measurements  were  as  follows:  Narrow  angle  radiometer 
measurements  yielding  q  were  made  viewing  across  the  exhaust  plume 
near  the  nozzle  exit.  (The  particle  flow  field  is  known  most  accurately  at 
this  location  in  the  exhaust.)  An  average  particle  temperature  along  the 
line  of  view  was  determined  as  described  below.  The  equation 

q=7acrbT4  (24) 


was  then  solved  for  P.  The  value  of  t  along  the  line  of  view  was  deter¬ 
mined  as  described  above.  The  effective  particle  emissivity  was  then 
determined  by  use  of  the  cylinder  curves  of  Figure  4. 


5 .  Specific  Values  of  Apparent  Emissivity 

With  the  value  of  r  determined  for  any  desired  view  through  the 
axis  of  the  particle  plume,  and  with  the  value  of  determined  as  de¬ 
scribed  above,  Figure  4  can  be  used  to  obtain  T ’  for  the  plume  segment 
encompassing  that  view. 


Each  segment  of  the  rocket  plume  is  treated  as  a  cylindrical 
particle  cloud  whose  axis  is  normal  to  the  view  line  from  the  surface 
element  through  the  center  of  that  segment.  The  value  of  t\  is  computed 
based  on  the  length  of  that  view  line,  f^,  per  Equation  (22).  This  value  of 
r.  along  with  c  are  used  with  the  curves  given  in  Figure  4  for  a  cylindrical 
cloud  to  obtain  a  value  of  €  a.  for  the  plume  segment.  The  value  of  ca^  so 
derived  is  taken  to  be  ca^  for  the  segment  surface  since  the  value  of  e  taken 
through  the  diameter  of  a  cloud  cylinder  represents  a  mean  value  for  that 
cylinder  }  * 

B.  EFFECTIVE  PLUME  TEMPERATURE 

The  particle  temperature  of  an  exhaust  plume  varies  with  particle  size 
and  particle  location.  Hence,  an  average  temperature  of  particles  along  a 
given  path  is  required  to  determine  the  heat  rate  from  these  particles. 

Since  heat  rates  measured  by  narrow  angle  radiometers  are  restricted  to 
a  small  field  of  view  of  a  cone  about  the  instrument  axis,  an  average  tem¬ 
perature  of  particles  along  this  axis  will  suffice  to  determine  the  heat  rate 
to  the  instrument.  However,  for  calorimeter  measurements,  the  field  of 
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view  is  hemispherical  (solid  angle  of  Ztt  steradians).  Therefore,  to  de¬ 
termine  a  heat  rate  to  a  calorimeter,  the  temperature  variations  of  the 
plume  in  lateral  as  well  as  longitudinal  directions  should  be  considered. 
The  longitudinal  variation  is  accounted  for  by  segmenting  the  plume  into  a 
number  of  slices.  However,  the  lateral  variation  is  accounted  for  by 
applying  an  averaging  correction  factor  to  the  average  temperature  of 
particles  along  a  pencil  view  through  the  centerline. 


Consequently,  determination  of  the  particle  plume  temperature  is 
reduced  to  computing  the  average  particle  temperature  along  a  given  line 
through  the  plume  axis.  Inasmuch  as  a  heat  rate  from  a  particle  cloud  to 
a  given  point  is  proportional  to  particle  temperature  to  the  fourth  power, 
to  particle  concentration,  particle  cross  section,  and  cloud  thickness,  and 
is  inversely  proportional  to  the  square  of  distance  between  the  point  of 
interest  and  an  individual  particle,  the  following  averaging  technique 
accounting  for  all  these  effects  has  been  used  to  determine  the  effective 
temperature, 


(25) 


The  correction  factor  rj  is  determined  by  taking  the  temperature 
variation  in  a  lateral  plane  of  view  and  averaging  it  over  all  view  lines 
contained  in  that  plane  weighted  by  lateral  variation  of  shape  factor. 

C.  TOTAL  RADIATION  FROM  PARTICLE  PLUME 

To  determine  a  total  heat  flux  from  the  particle  plume  to  a  given 
surface  element,  a  heat  flux  from  an  infinitesimal  slice  of  the  plume  must 
be  integrated  over  the  visible  region  of  the  plume.  The  integral  equation 
defining  the  total  heat  flux  is 

*  =  /  1  <a(r(v,KT4(Y)  (a?) dY  (26) 
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The  above  integration  may  be  performed  numerically  by  segmenting  the 
unshaded  portion  of  the  plume  with  a  number  of  path  planes  extending  from 
the  surface  element  to  the  plume,  as  shown  in  Figure  1 .  For  each  of 
these  segments,  the  mean  values  of  optical  thickness,  t.,  apparent 
emissivity,  T  ,  and  effective  temperature,  T. ,  can  be  determined.  Shape 

factors,  AF.,  from  the  surface  element  to  these  segments  can  also  be 

1  1 
determined  by  use  of  an  available  analytical  technique.  Then,  the  total 

heat  flux  from  the  plume  is  determined  as 


q  = 


T  (t.)T  AF. 
a.  i  i  i 

i 


(27) 


i 
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III.  ANALYTICAL  PREDICTIONS  AND  EXPERIMENTS 


A.  DESCRIPTION  OF  EXPERIMENTAL  TESTS 

1 .  Test  Models  and  Environment 

Experimental  data  were  taken  during  the  firing  of  two  aluminized 
propellant  rockets  under  simulated  pressure  altitude  conditions.  One 
rocket,  the  smaller  (A),  exhausted  through  a  contoured  nozzle  of  expansion 
ratio  19:1  which  had  an  exit  diameter  of  approximately  12  inches.  The 
second  rocket  (B)  exhausted  from  a  somewhat  similar  nozzle  of  expansion 
ratio  23.5:1  and  exit  diameter  about  26  inches.  The  chamber  temperatures 
and  the  percentages  of  aluminum  oxide  in  the  exhausts  were  representative 
of  those  for  typical  aluminized  propellants. 

The  firings  were  made  in  a  test  cell  arrangement  which  employed 
a  diffuser  and  an  auxiliary  pumping  system,  as  shown  schematically  in 
Figure  5.  Average  pressure  altitudes  of  110,000  feet  and  130,000  feet 
were  maintained  during  rocket  (A)  and  rocket  (B)  firings,  respectively. 
These  altitudes  were  sufficiently  high  to  permit  expansion  of  the  particle 
plume  into  as  large  a  plume  immediately  downstream  of  the  nozzle  exits 
as  would  occur  in  a  near  vacuum.  The  r ecirculatory  gas  flow  from  the 
boundary  of  the  exhaust  gas  plume  was  scavenged  by  the  auxiliary  pump, 
permitting  only  a  small  amount  of  convective  heating  in  the  area  at  and 
upstream  of  the  nozzle  exit. 

2.  Particle  Flow  Field 

A  particle  flow  field  was  developed  for  the  condition  of  each  of  the 
two  rockets  exhausting  into  a  near  vacuum  (Pc/Pe  =  10^)  by  the  computer 
program  of  Reference  2.  In  each  case  the  particle  size  distribution  was 
taken  to  be: 


Particle 

Number 

Radius, 

(microns) 

Percent  of 

ai2o3 

i 

0.79 

20 

2 

1.28 

20 

3 

1. 76 

20 

4 

2.44 

20 

S 

3.95 

20 

19 


This  distribution,  based  on  samples  collected  on  glass  slides  during 
simulated  firings  of  rockets  containing  similar  propellant,  results  in  a 
mean  particle  radius  of  2.04  microns.  Sample  flow  field  information  for 
the  two  extreme  particle  sizes  are  shown  for  rocket  (B)  in  Figures  6  and  7. 

3.  Instrumentation 

Instrumentation  for  measurement  of  thermal  radiation  was  of 
three  types:  (1)  narrow  angle  (5  degrees  solid  angle)  radiometers, 

(2)  black  surfaced  asymptotic  calorimeters,  and  (3)  spectrometers  which 
scanned  the  wave  length  range  0.3  to  4.5  microns.  Data  used  from  the 
rocket  (A)  firing  was  taken  by  narrow  angle  radiometers  which  viewed 
across  the  particle  plume  3  inches  downstream  of  the  nozzle  exit  and 
normal  to  the  plume  axial  centerline.  Data  taken  from  the  rocket  (B) 
firing  employed  all  three  types  of  instruments.  Figure  5  indicates  the 
instrument  locations  for  rocket  (B). 

B.  EXPERIMENTAL  RESULTS  AND  COMPARISON  WITH 

ANALYTICAL  PREDICTIONS _ 

1 .  Aluminum  Oxide  Particle  Emissivity 

The  aluminum  oxide  particle  emissivity  was  determined  as  de¬ 
scribed  in  Section  II  from  5  radiometer  measurements  from  the  rocket  (A) 
plume  and  from  7  similar  measurements  from  rocket  (B).  Values  are 
plotted  in  Figure  8  along  with  emissivity  measurements  made  on  bulk 
alumina^  at  elevated  temperatures.  The  reasonably  good  agreement 
between  the  two  suggests  that  the  particle  emissivity  approaches  that  of 
the  bulk  material  rather  than  the  low  values  of  sapphire.  This  seems  to 

corroborate  the  explanations  given  by  Lee  and  Kingery^  and  Kuby  and 
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Carlson.  A  value  of  particle  emissivity  was  also  determined  by  inte¬ 
grating  the  measured  spectral  intensities  over  wavelength  and  dividing  by 
the  integrated  black  body  intensity.  This  value,  =  0.3,  tends  to  con¬ 
firm  the  radiometer  measurements. 

It  is  seen  from  the  spectral  intensity  data  (Figure  9)  that  the 
plume  radiation  intensity  decreases  on  either  side  of  the  peak  wavelength 
at  a  rate  faster  than  that  of  black  (or  gray)  body.  Thus,  while  the  plume 
of  rocket  (B)  radiates  in  a  manner  somewhat  similar  to  a  gray  body,  it 
cannot  be  considered  a  gray  body. 
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Figure  6.  Particle  Concentration  Distribution 
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Figure  7.  Particle  Temperature  Distribution 
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Heat  Rates  to  Surface  Elements 


Using  the  method  developed  in  Section  II,  the  particle  flow  field 
for  rocket  (B),  and  a  particle  emissivity  of  0. 25,  predictions  were  made 
of  heat  rates  from  the  particle  plume  of  the  larger  rocket  to  the  calorim¬ 
eters  (surface  elements)  located  as  shown  in  Figure  5.  The  effect  of  the 
presence  of  the  diffuser  was  taken  into  account  approximately  in  the 
following  manner.  The  particle  concentrations  were  assumed  to  remain 
constant  downstream  of  the  point  at  which  the  theoretical  particle  limiting 
trajectories  intersect  the  diffuser  wall.  The  concentrations  were  taken  to 
be  equal  to  those  in  the  plane  of  that  intersection.  For  purposes  of  de¬ 
termining  t  and  T  for  those  segments  which  extend  downstream  along  the 
length  of  the  diffuser,  the  particle  plume  (constrained  by  the  diffuser)  was 
cut  off  at  three  nozzle  exit  diameters  downstream  of  the  nozzle  exit,  since 
downstream  of  this  point,  the  values  of  T  and  t  (and  therefore  <  )  would 
increase  such  a  small  amount  that  their  effect  on  the  contribution  to  q  to 
the  surface  element  from  that  segment  is  negligible. 

Table  1  compares  the  predictions  with  experimental  values  meas¬ 
ured  during  the  firing  of  rocket  (B).  The  agreement  is  considered  quite 
good  for  engineering  predictions  of  heat  rates. 


Table  1.  Comparison  of  Predictions  With  Experimental  Results 


=  o.i7 

q 
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The  relatively  large  differences  between  predictions  and  experi¬ 
ment  experienced  on  calorimeters  11  and  12  deserve  farther  comment.  It 
is  believed  that  convective  heating  due  to  recirculation  of  hot  gas  inside  the 
test  cell,  aluminum  oxide  slag  protruding  from  the  nozzle  wall  at  its  exit, 
and  reflection  from  the  supporting  bracketry  combined  to  increase  the 
measured  values.  These  effects  have  not  been  accounted  for  in  the  analyti¬ 
cal  predictions.  Including  these  two  measurements,  the  average  dif¬ 
ference  between  predictions  and  experimental  values  is  17  percent  of  the 
latter.  If  these  two  measurements  are  ignored,  this  value  is  reduced  to 
14.2  percent . 

C.  APPLICATION  OF  ANALYTICAL  METHOD 

It  is  the  opinion  of  the  authors  that  the  method  developed  here  can  be 
used  to  predict  radiant  heat  rates  incident  to  surfaces  located  in  the  region 
surrounding  the  plume  of  a  wide  range  of  particle  laden  rocket  exhausts. 
Examples  are  temperature  sensitive  spacecraft  surfaces  and  rocket  aft 
closures  and  components  located  near  the  exhaust.  Requirements  are 
that  the  particle  radiation  level  be  so  large  compared  to  that  of  the  gas 
that  the  latter  is  negligible  and  that  the  particle  material  and  size  dis¬ 
tribution  be  known  with  sufficient  accuracy  to  permit  development  of  the 
flow  field  by  means  of  a  two-phase  flow  computer  program.  The  particle 
emissivity  must  also  be  known. 

Predictions  should  be  suitably  accurate  for  design  purposes  except 
possibly  for  cases  in  which  the  surface  elements  are  oriented  such  that 
they  view  only  a  small  sector  along  the  outer  edges  of  the  exhaust  plume. 

An  example  of  these  surface  elements  is  the  outer  surface  of  the  nozzle 
wall  (Figure  1) . 

It  seems  clear  that,  in  this  instance,  the  approximation  of  a  plums 
segment  by  a  cylindrical  particle  cloud  is  not  valid.  Fortunately,  the 
heat  rates  from  the  plume  to  surface  elements  so  located  are  small  and 
are  usually  not  of  major  concern. 
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IV.  CONCLUSIONS 


An  approximate  method  has  been  developed  which  analytically  predicts 
radiant  heat  rates  incident  to  adjacent  surfaces  produced  by  a  rocket  exhaust 
containing  a  large  percentage  of  its  mass  flow  a_  s  aluminum  oxide  particles 
and  operating  in  a  rarefield  atmosphere.  Predi  ctions  can  be  made  with 
accuracies  suitable  for  engineering  design. 

The  emissivity  of  aluminum  oxide  particles  in  a  rocket  exhaust  appears 
to  be  between  0.  1  and  0.  3  and  is  not  inconsistert  twith  extrapolated  data 
for  bulk  alumina  at  a  temperature  in  the  vicinity^  of  3500°R. 

An  aluminum  oxide  particle  plume  radiates  in  a  manner  somewhat 
similar  to  a  gray  body.  For  wavelengths  on  eith  <er  side  of  that  at  which  the 
intensity  peaks,  the  rate  of  intensity  decay  is  considerably  greater  than 
that  of  a  gray  body. 
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